Soft swimming: Exploiting deformable interfaces for low-Reynolds number locomotion 
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Reciprocal movement cannot be used for locomotion at low-Reynolds number in an infinite fluid or near a 
rigid surface. Here we show that this limitation is relaxed for a body performing reciprocal motions near a 
deformable interface. Using physical arguments and scaling relationships, we show that the nonlinearities aris- 
ing from reciprocal flow-induced interfacial deformation rectify the periodic motion of the swimmer, leading to 
locomotion. Such a strategy can be used to move toward, away from, and parallel to any deformable interface as 
long as the length scales involved are smaller than intrinsic scales, which we identify. A macro-scale experiment 
of flapping motion near a free surface illustrates this new result. 
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Swimming microorganisms inhabit a world quite different 
from the one we experience. Their motion through the sur- 
rounding fluid occurs at very low Reynolds numbers (Re), a 
fact with two important consequences: (a) the only physical 
force available to produce thrust is drag; (b) since the fluid 
equations (Stokes equations) are linear and time-reversible, 
swimming motions symmetric with respect to time reversal 
(reciprocal motion) cannot be used for locomotion (the scallop 
theorem [1]). Biological swimmers, such as bacteria and sper- 
matozoa, overcome these two limitations by exploiting the 
anisotropic drag of slender filaments such as flagella and cilia, 
and actuating these filaments in a wave-like fashion [2, 3]. 

Motivated by the recent development of artificial swimmers 
[4], we pose here the following general question: are there any 
new low-Re swimming methods remaining to be discovered? 
In particular, can some simple reciprocal movements produce 
locomotion, thereby apparently violating the constraints of the 
scallop theorem? 

In this paper, we propose a new method for locomotion 
without inertia. Using theory and experiments, we show that a 
body performing a reciprocal movement is able to move near 
a soft interface: physically, the deformation of the interface 
provides the geometric nonlinearities necessary to escape the 
constraints of the scallop theorem. This strategy, which is 
relevant for sufficiently small systems, also implies that the 
time-reversible component of all swimmers (including biolog- 
ical organisms) can generate non-trivial flows and propulsive 
forces near soft interfaces. 

The physical picture for soft swimming arises from the 
asymmetries between different modes of motion near de- 
formable interfaces, and can be illustrated by considering the 
motion of a rigid sphere slowly translating below a free sur- 
face (Fig. 1). For motion perpendicular to the interface, the 
drag on the sphere is always increased by the presence of 
the interface; however, the magnitude and direction of the 
surface-induced component of the drag can be modulated by 
deforming the interface [5, 6]. When the sphere moves away 
from the interface (Fig. la), it displaces fluid with it, and the 
interface deforms toward the sphere, resulting in an increase 



of the surface-induced drag, Fdown > ^flat- In contrast, when 
the sphere moves toward the interface (Fig. lb), the free sur- 
face is displaced away from the sphere, and the component 
of the drag due to the interface is diminished, F up < Ffl at . 
The asymmetry, AF = Fd own — F up > 0> results in an av- 
eraged net force on the fluid, AF, away from the interface if 
the sphere is oscillated in a reciprocal fashion. (Alternatively, 
if the sphere is acted upon by an external force with zero av- 
erage, it will move toward the interface). For motion parallel 
to the interface (Figs, lc and d), surface deformation leads to 
a lift force on the sphere, perpendicular to its direction of mo- 
tion and directed away from the interface [7, 8], reminiscent of 
the lift force responsible for the migration of deformable bod- 
ies such as droplets [9], vesicles [10, 11] and polymers [12] 
away from boundaries. These effects are due to the nonlinear 
(geometric) coupling between the moving sphere and the de- 
formable interface, and they vanish when the interfaces do not 
deform. A similar coupling can be exploited for locomotion. 

Let us consider a body - a swimmer - that deforms in 
a prescribed and periodic fashion. On length scales larger 
than the typical size of the swimmer, the disturbance flow 
field can be modeled as a superposition of point flow distur- 
bances. Since a low-Re swimmer cannot exert any net force 
on the surrounding fluid, the instantaneous velocity field de- 
cays at most like a force-dipole (Fig. 2a); the swimmer is also 
torque-free, and therefore only symmetric force dipoles are 
allowed (stresslets [13, 14]). In three-dimensions, the dipole 
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FIG. 1: Breakdown of reversibility near soft interfaces: a sphere 
moving away from a soft surface experiences more drag than one 
moving toward the interface, Fdown > F up (a, b); a sphere translat- 
ing parallel to a soft surface experiences a lift force, F± , perpendic- 
ular to the direction of motion (c, d). 
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strength, p, is a 3-by-3 tensor with 6 independent components. 
The associated flow field decays as ~ 1/r 2 and is given by 
uP ~ s^r 3 ( — % + 3^r-) r, where r\ is the fluid viscosity, 
and r = x&i is the distance to the dipole. There exist two 
types of dipoles: parallel dipoles, in which the direction of 
the forces is aligned with that of the force gradient (Fig. 2b), 
and perpendicular dipoles, in which these directions are per- 
pendicular to one another (Fig. 2c). Furthermore, in general, 
the flow field can contain other types of singularities (higher 
moments of force distribution), with faster spatial decay, pre- 
dominantly the potential (source) dipole (Fig. 2d), of vector 
strength q, associated with a velocity field decaying as ~ 1/r 3 
and given by uf = ^_ (_£. . + 3 5i|i). 

In an infinite fluid, and by linearity of Stokes equations, 
the velocity of a swimmer whose far field is described by a 
superposition of stresslets and potential dipoles is given by 
u = M p + N q, where the mobility matrices M and N are 
functions of the instantaneous shape of the swimmer; here, we 
will take M and N as constants, a suitable approximation for 
bodies performing small- amplitude swimming motions. We 
also consider swimmers with sufficiently symmetric motion 
so that their orientation remains constant and ignore torque 
balance. 

We now examine the case of a swimmer located near an 
undeformable interface. The presence of the interface will re- 
sult in additional velocity for the swimmer, u w (h), where h 
is the distance to the interface, due to the image singularities 
required to enforce the correct boundary conditions at the in- 
terface [15]. The equation for the swimming velocity becomes 
u = M-p + N- q+u^/i), where the wall-induced velocity is 
given by u w (h) = 327r ^2 [(p xx +Pzz 2p yy )e y ^ (q x e x + 
q z e z + 2q y e y )] in the case of a flat free surface (illustrated 
by the dashed arrows in Figs. 2e-g); results for no-slip sur- 
faces are similar. Note that the magnitude of u w (h) scales 
linearly with the dipole strengths, p and q. In general, sur- 
faces of fixed shape will lead to an equation for the swimming 
kinematics given by ^ = A(r) -p(t) + B(r) -q(t). For recip- 
rocal deformations, we have p(t) = Po f(t), q(t) = qo/(£)> 
where / is an arbitrary periodic function of time: the swim- 
ming kinematics are described by an autonomous dynamical 
system, and no locomotion can occur on average, (u) = 0. 
This is the scallop theorem. 

We now consider interfaces that can deform due to the flow 
set up by the swimmer. Since wall effects depend on the 
shape of, and the distance to the interface, and since both can 
be modified by the swimmer's flow field, the relationship be- 
tween the local reciprocal actuation and the swimming kine- 
matics becomes nonlinear. 

Let us consider a free surface, with surface tension 7, and 
show that motion toward, away from, and parallel to the sur- 
face can be devised. For motion toward the surface, we only 
need to consider parallel force-dipoles aligned with the sur- 
face, p = p(t)e x e x (Figs. 2h-i). We assume small Capillary 
number to capture the first-order effect of surface deformation 
(p/jh 2 <C 1). The normal deformation of the interface (mag- 
nitude S) is set by the balance between capillary and viscous 
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FIG. 2: Physics of soft swimming, (a) The lowest-order flow around 
a low-Re swimmer decays as a force-dipole; (b): parallel force- 
dipole (p xx > 0); (c): perpendicular force-dipole (p xy < 0); (d): 
potential (source) dipole (q x > 0); (e): p xx > is attracted by a flat 
surface; (f): p yy > is repelled by a flat surface; (g): q x < moves 
parallel to a flat surface; (h): p xx > deforms a soft surface in a 
manner that enhances the attraction; (i): p xx < deforms the same 
surface in a manner that reduces the repulsion; (j): a reciprocal p yy 
is attracted by a soft surface; (k): a reciprocal p xy is repelled by a 
soft surface; (1): a reciprocal combination of p xy and q x moves away 
from, and parallel to a deformable surface. 



forces at the interface, jS/h 2 ~ —p/h 3 , hence S ~ —p/jh; 
note that the surface deformation can be calculated analyti- 
cally for small Capillary numbers, as plotted in Figs. 2h-l. The 
wall-induced velocity on the dipole, due to the image system, 
now has a magnitude u w = e y ■ u w ~ p/rj(h + S) 2 , and can 
be linearized around ho (the starting position of the swimmer) 
to obtain u w ~ p/rj(h ) 2 -p5/rjh% ~ p/rj(h ) 2 + p 2 /^vK- 
The swimming kinematics in the y direction are then given by 
dh/dt ~ Mp(t) + p(t) 2 /jrjhQ —p(t)/rjhQ. For reciprocal 
swimming p(t) = po cos out, we get a net locomotion toward 
the surface, with average swimming velocity (u y ) ~ Po/lV^o 
[19]. In practice, the dipole arises from the unsteady actuation 
of a swimmer of size L with velocity U ~ coL , and therefore 
po ~ r]L 2 U, so that (u y )/U - (rjU/j)(L/ho) 4 : for suffi- 
ciently large deformations 77/7/7 ~ 1 and confined swimmers 
L/ho ~ 1, we obtain a net speed on the order of the actuation 
speed (u y ) ~U [20]. 

The physical interpretation for this motion toward the sur- 
face is illustrated in Figs. 2h-i. When p > 0, the swimmer de- 
flects the surface toward it, enhancing the wall-induced attrac- 
tion (Fig. 2h, dashed arrow) by a small amount (empty arrow); 
when p < 0, the surface is deflected away from the swimmer, 
reducing the wall induced repulsion (Fig. 2i, dashed arrow), 
by a small amount (empty arrow). Both deformation-induced 
effects are directed toward the surface, and therefore do not 
average to zero for a reciprocal actuation. The results are sim- 
ilar for a parallel force dipole in the direction perpendicular to 
the surface, p = p(t)e y e y (Fig. 2j). 

Locomotion away from the surface can be obtained by 
considering perpendicular force dipoles, of the form p = 
p(t)(e x e y + e^e^) (Fig. 2k). In that case, the surface de- 
formation is asymmetric, and it induces a lift force which 
does not average to zero with p [5]. The scaling for the 
time-averaged swimming speed is similar to the previous sec- 
tion, except that it is now directed away from the free surface, 
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FIG. 3: Examples where soft swimming can occur: (a) interface be- 
tween two liquids (or between a liquid and a gas), (b) elastic mem- 
brane separating two liquids, and (c) interface between a liquid and 
an elastic gel (see text for notations). 



(u y ) ~ -Poftvht 

Moving parallel to the deformable surface is more subtle. 
Indeed, since the first effect of surface deformation is always 
directed perpendicular to the undisturbed surface [5], the only 
way to move parallel to the surface is to use motion toward 
or away from the surface to rectify an otherwise time-periodic 
and spatially dependent parallel motion (such as the one illus- 
trated in Fig. 2g). This can be achieved, for example, by com- 
bining a perpendicular force dipole, p = p(t)(e x e y + e^e^), 
with a potential dipole, q = q(t)e x (Fig. 21). Since the flow 
field due to the potential dipole decays faster (by one order 
of magnitude) than that due to the force-dipole, the first ef- 
fect of surface deformation will be the establishment of the 
lift force due to p. The swimming kinematics are therefore 
given by Ah /At ~ —p(t) 2 /^r]h A in the direction perpendic- 
ular to the surface, and Ax /At ~ —q(t)/r]h 3 in the parallel 
direction. Writing h(t) = ho + h(t), with h(t) <C ho, we 
have Ax /At ~ —q/v^o + qh/r)h\. For a reciprocal motion 
with p(t) = po cos(ut + tp) and q(t) = qo cos(ut + cp), time 
integration leads to h(t) = —p 2 ) (t+sin(2u;t-\-2(p)/2u;)/'yr]hQ, 
and therefore the average swimming speed in the x direction 
is given by (u x ) ~ — p^qo sin (p/uj^rfho- Thus swimming 
occurs both away from the surface and in the + or — x direc- 
tion, depending on the sign of sin cp. Similarly, combining a 
force dipole in the e x e x or e y e y direction with a source dipole 
along e x leads to combined motion toward and parallel to the 
surface. 

The examples above demonstrate that a free surface can be 
exploited to obtain locomotion in all directions using recipro- 
cal forcing. To analyze the effectiveness of these swimming 
strategies, three types of biologically-relevant soft surfaces 
can be considered (Fig. 3): the interfaces between two liquids 
(Fig. 3a), an elastic membrane separating two liquids (Fig. 3b) 
and the interface between a liquid and an elastic gel (Fig. 3c). 

Consider a fluid with viscosity rj and density p, in which all 
length scales L (swimmer size, a, distance to the surface, h) 
are similar. For soft swimming to occur, the local reciprocal 
actuation of the swimmer has to generate flow speeds, U, suf- 
ficiently small for the Reynolds number, Re = pUL/rj, to be 
small, yet large enough to induce significant deformation of 
the interface, i.e. a large (a) Capillary number, Ca = rjU/j, 
for the interface between two fluids with similar viscosities 
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FIG. 4: Experiment to generate a net flow by reciprocal flapping mo- 
tion below a free surface; (a) principle of the experiment; (b) velocity 
and vorticity fields over one flapping cycle (red is positive vorticity); 
(c) dimensionless kinetic energy over one cycle as a function of the 
flapping Capillary number, Ca (insert, photograph of the free surface 
deformation). 



(rj, rj*) and densities (p, p*), (b) dimensionless membrane- 
viscous number, Me = t]U/Ka, for two fluids separated by 
an elastic membrane (thickness t, area modulus Ka, den- 
sity p m , bending stiffness Kb), (c) elasto-viscous number, 
El = r]U/LE, for an elastic gel (Young's modulus E, Pois- 
son's ratio, v g , density p g ). By combining Re < 1 with (a) 
Ca > 1 , (b) Me > 1, (c) El > 1, we find that soft swimming 
is expected to be effective for length scales below (a) L < if, 
where if = rj 2 /fry is an intrinsic visco-capillary length scale 
(the Ohnesorge length), (b) L < i m , where i m = t] 2 /pKa 
is the intrinsic membrane- viscous length scale, (c) L < i g , 
where i g = ^/r} 2 / pE is the intrinsic elasto-viscous length 
scale [21]. 

In order to probe the biological relevance of these ideas, 
we consider an environment where motility is crucial, namely 
the female reproductive tract in humans [16]. In this case, the 
fluids involved (e.g. cervical mucus) have viscosities in the 
range 77 ~ 0.1 — 10 Pa.s [17], leading to values of if between 
fractions of millimeters to meters. A similar range for i m is 
obtained in the case where such fluids are separated by a lipid 
bilayer [18], and for i g when the fluids are near an elastic gel 
such as a cross-linked actin network [18]. In all cases, the 
intrinsic length scales are large, and therefore biological flows 
arising from time-reversible motion could be substantial. 

To further illustrate the breaking of the scallop theorem near 
soft surfaces, we exploit the increase of the intrinsic scale, if, 
with the viscosity, and perform a macro-scale experiment near 
a free surface. Instead of a swimmer, we consider a flapper 
fixed in place, for which the problem of reciprocal swimming 
is replaced by the dual problem of reciprocal pumping, i.e. 
creating a non-zero flow from a time-periodic forcing. 

A thin plate, the flapper (size 2 cm x 7 cm), was mounted 
to the shaft of a servo motor and rotated over a range of 90° 
(45° on each side) at flapping periods of 3-200 s (Fig. 4a). 
The flapper was immersed in a high-viscosity silicone oil 
(77 = 60 Pa-s, if w 50 m) with the tip of the flapper positioned 
2 cm below the free surface; alternatively, the free surface 
could be covered by a solid acrylic sheet to prevent surface 
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deformation. A laser sheet illuminated poly amide beads (di- 
ameter, 93 pm) dispersed in the oil, and images of the beads 
were captured after every flapping cycle using a digital SLR 
camera with 3072 x 2048 pixels per image. The displacement 
field per flapping cycle was calculated from sequential images 
using the MatPIV toolbox for Matlab. 

The experimental results are presented in Fig. 4b-c. When 
the surface is free to deform, the reciprocal flapping motion 
induces a net flow. An example of the vorticity contours (red 
is positive vorticity) and flow field over one flapping cycle is 
displayed in Fig. 4b: the net flow is induced along the flapper, 
away from the free surface; when the free surface is replaced 
by a solid lid, this flow disappears. Measurements of the ki- 
netic energy in the flow (integrated over the whole system) are 
shown in Fig. 4c (nondimensionalized by flapper length and 
flapping frequency). When the surface is solid (hollow sym- 
bols), no flow occurs on average, whereas when the surface 
is free to deform, the average flow possesses finite kinetic en- 
ergy (filled symbols), demonstrating that reciprocal actuation 
near a soft surface can be used to escape the constraints of the 
scallop theorem. 

In summary, we have shown that soft interfaces allow 
propulsion driven by reciprocal motion. This new strategy for 
low-Reynolds number swimming means not only that primi- 
tive swimming methods become effective near a soft interface, 
but more generally that the usually-neglected time-reversible 
component of all swimmers (including biological organisms 
and appendages) can generate non-trivial flows and propul- 
sive forces in such settings. The symmetry-breaking arises 
from the normal stress boundary condition: changing the sign 
of the forcing by a reciprocal swimmer changes the sign of the 
interfacial normal stresses, thereby changing the interface cur- 
vature, and its shape. The strategy is relevant for sufficiently 
small systems. Moreover, for swimmers of size L at a dis- 
tance h from a surface, the swimming efficiency is expected 
to scale as £ ~ (L/h) n (n > 1) and will therefore be most ef- 
ficient for confined geometries where L ~ h, for example for 
the motion of swimmers in thin films and membranes, or near 
surfaces. Similar results are expected for the dual problem of 
generating feeding currents from reciprocal forcing. Possible 
extensions of this work include investigating surface-induced 
reorientation, rotational swimming, and a closer examination 
of cases in which surface deformations are large. Together 
with the current work, this would provide a comprehensive 



framework for understanding the consequence of soft inter- 
faces on locomotion and transport on small scales. 
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